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Abstract 



We consider a wide class of two-dimensional models as gauge the- 
k^ I ories, Gross-Neveu model, 0{N) and C-P^~^-like models using a for- 

;h ' malism based on the introduction of bilocal fields that permits to 



perform easily the large-N expansion of this set of models in a unified 
and general way. We mainly discuss the SU{N) gauge field theory 
minimally coupled to fermionic plus bosonic matter in the fundamen- 
tal representation, and we obtain within the path integral approach 
exact equations for the particle spectrum, also in presence of renormal- 
izable polynomial potentials. Finally, we discuss the correspondence 
between this new approach and the one previously used in the context 
of the 0{N) vector models. 



E-mail: cavicchi@bo.infn.it / cavicchi@nbivax.nbi.dk 



1 Introduction 

One of the most interesting and still not fully understood problem of Quan- 
tum Field Theory concerns the calculation of the mass of relativistic bound 
states. This problem is even more complicated to solve in the case of QCD, 
the SU{3) gauge theory of quarks and gluons, that describes the strong in- 
teraction physics where the basic constituents are confined and perturbation 
theory cannot be used. The expectation is that the low energy spectrum 
consists of colorless mesons and baryons, but up to now the only available 
method seems to be based on numerical calculation using the lattice theory 
formulation of the theory. 

The large-A^ expansion technique proposed by 't Hooft |l[] several years 
ago seems to be the most promising approach to obtain analytic results on 
the hadron spectrum. It is well known that in the limit in which the num- 
ber N of colours becomes very large the theory becomes much simpler, in 
the sense that only planar Feynman graphs survive. Furthermore, when 
N ^ oo the theory only contains colorless, stable and noninteracting mesons 
with two-body decay and scattering amplitudes proportional to -4^ and -^, 
respectively (there is also a way of describing baryons as solitons of the ef- 
fective Lagrangian in the large N limit, but this will not interest us in this 
context). 

Unfortunately, all efforts to try to solve four-dimensional large- A^ QCD 
have failed up to now. The main problem is to find the semi-classical con- 
figuration, namely the master field Q, from which the action is dominated 
in the large- A^ limit and whose fluctuations around the vacuum should give 
the particle spectrum of the theory. 

The main reason of this failure is due to the fact that no method has 
yet been found to solve matrix models for a space time dimension D > 2. 
It is well known that the large N expansion can be explicitly performed in 
vector models as for instance the 0{N) vector model (see for instance ref. 
and references therein), the two-dimensional CP^~^ model and also in 
QCD2 @], [§ with matter in the fundamental representation of the gauge 
group where the gluon field can be eliminated by using its equation of motion. 

In particular this set of models are solved by means of two slightly dif- 
ferent methods. The vector-like models are solved in the large N expansion 
by introducing a local composite field and by explicitly integrating over the 
fundamental fields (See for instance Refs.[Q and |^). The QCD2-like mod- 



els are characterized by the fact that the number of components of the non 
abehan gauge field goes to infinity in the large N expansion. The gauge field 
in two dimensions has no physical degrees of freedom and therefore can be 
eliminated by using its classical equation of motion. In this way one gets a 
non local Coulomb interaction that is quartic in terms of the vector-like mat- 
ter fields. The theory can then be solved in the large N limit by introducing 



a bilocal composite field as discussed in Refs. and |23 . 

In this paper we consider a general two dimensional gauge theory with 
matter transforming as the fundamental representation of the gauge group 
(vector-like matter) and we solve it by introducing bilocal colorless fields in 
order to explicitly perform the large-iV expansion. The introduction of a 
bilocal field is essential because the gluon field has a number of components 
that goes to oo in the large N expansion. We also show how our formalism 
reduces to the usual formalism of pure vector like models in the limit where 
the gauge coupling constant g^ -^ 0. 

The work is organized as follows. In section 2 we write the Lagrangian for 
a SU{N) gauge theory minimally coupled to fermionic and bosonic matter 
in fundamental representation in the light-cone gauge. 

In section 3 we consider the theory with only scalar fields present, at 
first in the free massive case, then adding with a self-interaction potential 
t/(0V)j we solve it in the large N expansion and we show how to obtain from 
our solution the known results for the 0{N) model in the case in which the 
gauge coupling g is set to 0. In the last subsection we constrain a (massless) 
interacting scalar field to satisfy 0|(x)0j(x) = N/2f, to find some similarities 
with the CP^"^ model. 

In section 4 we consider the theory with fermionic fields, already studied 
in |2^, but adding current-current couplings; this kind of model is also 
compared with the extension of the CP'^~^ model with quarks. 

Finally in section 5 we study the full gauge theory, (fermionic plus bosonic 
matter in fundamental representation of the gauge group) showing how to 
obtain the corresponding of the 't Hooft equation p| for all particle sec- 
tors (b-b, b-f and f-f bound states) and examining the effect of adding any 
renormalizable and polynomial self-interaction term. We write exact integral 
equations, that reduce themselves, in the case of zero gauge coupling, to the 
equations very similar to those previously found for the spectrum of the mas- 
sive Thirring model; in this case we find the exact spectrum of states. We 
also examine the case of nonzero Yukawa couplings between the fermionic 



and the bosonic SU{N) fields. 

2 Two-dimensional gauge theory with mini- 
mal coupling 



We start with the Lagrangian 



C 



-f;,f'^^^ + ^(^ 



mo)V + pM0rP''0)-"^o0V 



(2.1) 



where F, 



flU 



d,A, 



All matter fields are in the fundamental representation of the gauge group 
SU{N) and we have chosen T"" to be a Hermitian matrix. 

Then we move to the light-cone coordinate system Q, and we interpret 
the x~^ coordinate as our 'time' coordinate. The Lagrangian is particularly 
simple in the light cone gauge A°i = 0, because the nonabelian part of the 
field strength and his quadratic coupling with the boson field vanish: 



jCu, = ^id^Alf + i,{%f - mo)^ + 0t(-92 - m, 



g^Air 



(2.2) 



with 



r_ = ipT^-f.iP + icP^T"" d. 



(2.3) 



The equation of motion for the field A°^ does not contain any 'time' derivative, 
so we can eliminate it to get an effective (non local) Action involving only 
the matter fields. 



S, 



eff 



d^xcPy I il){i(p — mQ)'ip + 0^(— 9^ — rrig 



3 Bosonic theory 



^xy + 2^0 /}2 '-'- 



(2.4) 



*The conventions are the same as in fE3| 



3.1 The 't Hooft equation for the bosonic sector 

We now want to study the theory in the case in which only scalar fields are 
present; the case with only fermions in the theory has already been treated [^ 
with the same methods. The model without field self-interaction, some- 



time also called "Scalar QCD", has been extensively studied ||14[ using the 
same techniques developed by 't Hooft; we now want to derive the same re- 
sults by our new approach. We rewrite the interaction term in the following 
form: 

S^nt = \gljdPxdPy (</)t(a;)T" C </)(x)) G{x - y) {(^\y)T'^ cf- 0(y)) 

(3.1) 
with the Green's function G defined by 

- dtG{x - y) = 6^^\x - y) (3.2) 

so that in momentum space G{k) = 1/k'^. In the following, for the sake 
of simplicity in the notation, we will use di and 82 instead of d^- and dy~ , 
respectively. 

By using the relation, that is valid only for fields transforming according 
to the fundamental representation p3[ 



zZ'^AB^CD — ^BC^AD — -TZ^AB^CD (3.3) 



we see that it is convenient to introduce as in [EH] a bilocal colour singlet 
field 

A 

The 0(1/N) term in eq. (|3.3|) is completely irrelevant to our purposes, and 
we will neglect it in the following. 

With the help of ( p.3|) we can rewrite ( |3.1|) : 



■^^ ■^— > 



C,nt = \9lG.,ya^Hd^id2ay^ (3.5) 



The jacobian of the change of variable from (jrcj) to a can easily be computed 
in the large-N limit (see [^, for example), and gives a contribution to the 
action that can be written formally as follows 

5S[a] = -iNTr\oga (3.6) 



Now, rescaling a — >■ Na and sending N ^ oo keeping g^N = g^ fixed, we 
see that the action is dominated by a very simple saddle point. The action 
in the large-N limit is then given by: 



S/N = -iTr log a + / d^xd^y 
in configuration space and by 



(3.7) 



S/N = -iTr\oga+ f 



d'^p d'^q 

(27r)2 {27iy 



ip'-ml)aip,q)i2nf5^'\p + q) 



g^ f d^pi d^p2 d^k (2pi + A;)_(2p2-A;)- 



+T 



(27r)2 (27r)2 (27r)^ 



fc? 



(y{pi,P2)(y{k - P2, -k - pi] 

(3.8; 



in momentum space.Q 

The master field ctq is the bilocal field configuration that dominates the 
functional integral for large A^. It satisfies the saddle point of the action 

= 0; As explained in Ref. ||2^ translational invariance imposes to the 



<5cr 



(T=cro 



^ We must give a suitable definition of —d_ or, in momentum space, of 1/k^. We can 



easily see that the equation 



^dlG{x ^ y) = S'-^\x - y) 



has the general solution 



G{x) = S{x+) 



-2l^ 



Bx- 



A 



The terms A and Bx^ are the zero- modes of the operator d^. As it has been pointed out 
by Einhorn |^, the B term corresponds to a background colored electric field, and without 
loss of generality it can be taken as zero; furthermore, the constant A can be gauged away, 
so that we can conclude that it contains no physics. In his work, 't Hooft 1^] regularized 
his integrals by cutting off all momenta \p- \ < A where A was a constant to send to zero 
at the end of all calculations. It can be demonstrated that this procedure corresponds to 
the choice for the constant A to be equal to I/ttA so that in momentum space 
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1 

"2 


1 1 


kl 


_(fc_+2e)2 (fc_-ze)2_ 



In the following, for those "historical" reasons, we will keep this regularization prescription. 



master field to be local i.e. 

Mp,<i) = Mp)i'^^ys'^'\p+<i) (3.9) 

A simple algebra permits to rewrite the master field equation as follows 

(^p' - ml - g' J -0^^^^^^^ao{p + k)yo{p) = t (3.10) 

Defining 

5^' = //^^o(A;) (3.11) 

and using the following ansatz 

Mp)=KP^-m^ + '^ip)+^^y^ (3.12) 

with m'^ = ml + Sm^, we get the equation for T{p) 

J.. . 2 ( d'k (2p_ + k_)' - e i 



(27r)2 e_ {p + A;)2 - m^ + r(p + fc) + ie 

that has the solution |0| 

r(p) = £! _ £!k^ (3.14) 

In the case of scalar fields, unlike the fermionic case, we needed to renor- 
malize the mass (the self-mass is affected by the ultraviolet divergence of the 
boson loop). 

The particle spectrum of the theory is obtaining by studying small fluctu- 
ation around the master fleld. Putting a = ctq + ~m^^i ^^ ^^d the effective 
action for the fluctuations by keeping only the quadratic part in (5o", from 
which we get the equation for the fluctuations: 

. / N • 2 / X / . [ (Pk (2p + A;)_(2g- A;)_ ^ ,, , . 

5a(p, q) = ig ao{p)ao{-q) J -^^ p da{k + p, -k + q) 

(3.15) 
The bound state equation can now be obtained by simply following the same 
procedure as in Refs. 0, |2^. One deflnes 5a{r, s) = 6a{p + q, {p — q)/2) and 
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after the integration of both sides of the previous equation over ds +/27t, the 
previous equation becomes (let us choose r_ > 0): 





r M2 M2 g^ 1 
1 \ r 


V9(s_) 






_2s_ + ^| 2|s_-^ ttA ^ 




= 


/ /•*-+^ rfA;_ (2s_+r_ + A;_)(-2s_+r_ 


-A:_) 


27rA_-l^ A;2 2|s_ + ^|2|s_ - ^ 




The (p field is simply defined as 






V{s_)=l^j~<r.s) 





V9(s_ + kJ^ 



(3.16) 



(3.17) 



(the r-dependence has been left implicit) and by definition r is the total 
momentum of the two-parton system. 

The A-dependent part cancels applying our regularization 



/ -^ <^(^- + ^-) = jV{s-) + P j ^ ^^'- + ^ 



(3.18) 



Rescaling and shifting the variables, we then easily get the 't Hooft equation 
for the bosonic sector (r^ = 2r_r i ): 



r'^ip{x) 



M2 M2 

X {1 — x) 



9" r^ f^ dy {x + y)i2 ~ X - y) 
2x2(1 -x) 



TT Jo [y — xy 



v{y) 

(3.19) 



3.2 Self-interacting scalar field 

In this subsection, we want to analyze the effects of adding a potential 
V{(j)^(j)). \n D = 2 dimensions we can choose an arbitrary potential: 



M 



Ck 






(3.20) 



without losing the renormalizability of the theory. 



After the introduction of the bilocal field a potential term adds to the 
action a local term of the form 






(3.21) 



After the rescaling a —>■ Na we see that in the large A^ expansion we have to 
keep constant the quantities Ck = CkN^~^. For simplicity we will only study 
a Y^'P^^py potential, and let it be f = C2 = C2N. 
The new equation for the saddle-point is quite simple: 



222 
P -^Q-g 



(fk {2p_ + k_ 



-c7o{p + k)-v J 7^^^o(^) 1 Mp) 



(3.22) 

and it can be seen that the effect of the potential term on ( p.lO| ) is just to 
change the definition of the mass shift Sm^: 



6m^ = {g^ + v) I 



(Pk 

(2^ 



ao(A;) 



(3.23) 



We can then easily compute, as in the previous section, the equation 
for the fluctuations around the master field, getting the equation previously 



written by Ambj0rn Q (M^ = 



m 



gV^) 



r'^Lp{x) 



M2 



X 



+ 



M2 



X) 



i^{x) 



-p 



+ 



V 



dy {x + y){2-x-y) 



(y-x)2 
1 n 



v{y) 



4:71 X(l — X) 



-p 



2x2(1 -x) 
dy^iy) (3.24) 



that is equal to eq. ( 3.19|) with the addition of the last term coming from 
the quartic potential. 

If ^f = the gluon field, that is a matrix, decouples from the matter and 
we are left only with vector-like fields and it is well known that one needs 
to introduce only a local composite field (see for instance Ref. [|]) and not a 
bilocal one. Therefore if we set g = the two methods must agree and this 
is what we are going to check in the following. 

Let us start from the master field equation in eq. (|3.10| ) for the case 
g = whose solution is given by the expression in eq. ( |3.12| ) with T{p) = 0. 



Integrating both terms of eq. ( p.l2| ) over p and defining 



— Ioq = 5m 



(27r)^ 



rf^o(^) 



(3.25) 



we get, after a Wick rotation, the following equation 

iv J (27r)^ —p^ + m? — iao 



(3.26) 



that is identical to equation (2.25) of Ref. with 4/ = v, Kf = Kn = 0. The 



r.h.s of eq. ( |3.26|) is divergent and must be regularized. The renormalized 
equation is obtained by extracting the divergent piece from ctq ( this divergent 
is due to the fact that (Tq is the vacuum expectation value of a composite 
local field where the two constituent fields are taken at the same point) and 
having it to be cancelled by the divergence appearing in the r.h.s. of eq. 
( |3.26| ). After this procedure one obtains a renormalized gap equation ( See 
eq. (3.15) of Ref. 0). In conclusion we have shown that, ii g = the master 
field equation reduces to the gap equation of the vector models. 
Let us consider now the fluctuation equation for g = 



m 



x(l — x) 



^^"^=i4T^''r^'^^'^ 



(3.27) 



whose solution is given by 



(p{x) 



V 



An r^ x(l — x) 



m- 



:P J dyip{y) 



(3.28) 



By integrating over x we get the consistency condition 



dx 



— P 

4:71 Jo r^xfl — x) 



m 
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2 ^i^^^K^') 



(3.29) 



that gives three cases, depending on the region in which we are looking for: 
i) r^ > Am?: (B.29I) becomes 



47r^2 



4m2 2 



1 1 + a/1 
log- 



im'^ 



1-a/1 



4m,2 



(3.30) 



or, parametrizing r = 2m cosh 7, 7 > : 

'^^^^ = -^ (3.31) 

27 Aum^ 

that has an unique solution iff f > Anm"^, where we have a resonance of mass 
^2 = f2 ^j^]^ ^2 ^^g^^ satisfies (p.30| ). 



ii) < r^ < 4m^: ( p.29|) becomes 

i; 4 



47r^2 /4mi_^ ./4l#-l 



Atg^=^= (3.32) 



^2 



that is, parametrizing r = 2m sin 6*, < < 7r/2 : 

sin 26 V 

29 ~ ~ Atxw? 



(3.33) 



This condition can be satisfied only if —1 < v/Airm'^ < (see Fig.l), and we 
have a bound state determined by ( p.32|) . 

iii) r^ < 0. Obviously, if we find tachyonic solutions of ( |3.29| ) the theory 
is not defined. Writing r^ = — 4m^z/^, i/ > we get 



1 



Airm'^ y/T 



V 



log(z/ + Vl + 1/2) (3.34) 



that admits solution iff v/ATrm? < — 1 (see Fig.2), so that this range of values 
for the coupling constant must be avoidecfj. 

We can find a correspondence between the approach proposed here and 
the methods that people used in the past. In the usual procedure for the large 
N expansion in the context of the 0{N) vector model one gets a quadratic 
term for the fluctuations around the solution of the gap equation that can 
be written as (eq. (2.29) of Ref. g ) 

si% = lJ-0^A^nk'M-k) (3.35) 

r(fc^) (not to be confused with our r(p_)!) is defined by eq. (2.30) of 
Ref. 0. The equation of motion for the Fourier transform of the fluctuation 

^ This result agrees with the spectrum obtained using semiclassical methods by Abbott 
pOl , who pointed out that the condition v/ATrm^ > — 1 was necessary in order for a vacuum 
of the theory to exist (see also [O for a more detailed discussion) . 
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field a{x) gives r(/i;^) = (with r(A;), in the D = 2 case, that is given by eq. 
(3.30) of Ref. 0). This condition is just equal to eq. (p.30|) , Wick rotated 
to Euchdean space and with v = 4/. 

Another more direct way to find the connection between our approach 
and the one in terms of a local composite is the following. If we consider the 
equation of motion for the field 6a obtained from the action (|3.8| ) plus the 
self-interaction term, the Fourier transform of ( 3.21|) , setting g = we have 



r (Pk 
5a{p, q) = -ivao{p)ao{-q) J 7^:^^(^i^ + p, -k + q) (3.36) 

It can be seen from (|3.36| ) that the field (7o^{p)6a{p,q)aQ^{—q) depends 



only on p + g and not by p and q separately. Therefore, we can identify the 
local field a of Ref. with this combination of fields: 

a{p + q)= a^^{p)6a{p, q)(jQ^{-q) (3.37) 



Substituting ( |3.37| ) into ( p.36| ) we can eliminate a and we immediately 



get eq. (|3.29| ) (after the introduction of a Feynman parameter and the inte- 



gration over k). 

The above discussion implies that the bilocal approach that is essential to 
treat large N theories with matter transforming according to the fundamental 
representation of the gauge group ( vector-like matter) reduces in the limit 
of zero gauge coupling constant to the standard large N approach used in 
the pure vector models. 

Finally it can be easily seen that a complete potential of the form ( p.20| ) 
has just the effect to change the equation for the renormalization of mass 



//7g.-o(P) + E7^(/|^^oW)" (3.38) 



Here we have to define a new renormalized scalar self-coupling to be used in 
the equations for the fluctuations 

.„ . V" (/.„) ^ I j^ (/|f5-o(P))'" (3^39) 
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This procedure is in agreement with the one presented in Ref. 0. 

From the above discussion, we have seen that the bilocal field formahsm 
is a more general method than the one based on a local composite, that can 
be applied both to the pure vector models and to models containing a gauge 
field whose number of components goes to oo as A^ — ;► oo. 

3.3 A CP^^ -like model in bilocal field approach 

The model that we will discuss in this subsection is a massless scalar field, in 
the fundamental representation of the SU{N) gauge group, minimally cou- 
pled to a gauge field in adjoint representation and subject to the constraint 
0|(a;)0j(a;) = N/2f. The main difference with respect to the CP^~^ model 
is the presence of the kinetic term for the gauge field and the fact that the 
gauge field is not abelian, so that its number of components grows with A^. 
We therefore start with the model 

C = --^F^^F'^^^ + (Z}^</))t(Z}^0) (3.40) 

The potential for the scalar field is replaced by the constraint 

0l(x)0,(x) = iV/2/o (3.41) 

where /o is a dimensionless constant. We can express the constraint ( p.41| ) by 
introducing a Lagrange multiplier in the action, so that the partition function 
will take the form 

Z= /'l)aPAe*^M-*/'^'^^-('^^^-^/2/o) (3 42) 

Where now A^. is a (local-field) Lagrange multiplier not to be confused with 
the constant previously used in the definition of 1/fc^. After the usual rescal- 
ing 0" -^ Na we can obtain an action 5* [A, a] that can be solved in the large- A^ 
limit via a saddle point approximation. 
The saddle point is obtained by imposing the two conditions: 

^ = (3.43) 

OCT cr=(TQ,A=AQ 

f = (3.44) 

OA A=Ao,cr=cro 
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They imply the following equations: the one for a 

(3.45) 
and 

"^'^^aolpi,/ -pi) = -V (27r)25(2)(/) (3.46) 



^a^\p,q)=p\27^f5^'\p + q)-\o{p + q)-g' j ^^^y ^2 MP + k) 



The main difference between (|3.43|) and ( p. 441) is that A, being a /oca/ field, 



has a constant vacuum expectation value. So we can write Ao(/) = A ■ {2nY5^'^' (/) 
and the previous equations become 



p. ^ 2 1 d^k (2p_ + k^f - e i 



and 



or 



(27r)2 kl {p + kY-m? + T{p + k)+ie 

1 /• rf2p 



2/0 7 (27r)2p2_^2 + r(p)+^g 



(3.48) 



with the usual ansatz ( |3.12| ) for ctq) and with m^ = A + (5m^. We needed a 
renormalization of / (27r//o = 27r// + log(A2//ig)) and /ig is the renormaliza- 
tion scale. Eq. ( p. 49] ) is expression of the asymptotic freedom of the model, 
because the coupling constant goes to zero as the scale /Xg goes to infinity. 
Comparing ( p.48| ) with (40) of [0, one gets some similarities; in fact, the 
only difference is the replacement of the mass parameter m"^ with m"^ — g'^/ir: 
even in this gauge-interacting (and then nonlocal) version, the dependence 
on the cutoff of the coupling / remains the same, and so it happens for his 
'asymptotic freedom' property. 

We then can see how the model with a potential is related to the one with a 
constraint. In both cases we get the same equation ( compare eq. ( p.4?| ) with 
eq. (|3.13|) ), but in one case the bare squared mass ml enters directly in the 
Lagrangian while in the other case is generated by the v.e.v. of the field A^;. 
When the kinetic term for the gauge field is absent ( as in the real CP^~^ 
model ) we simply have the gap-equation ( p.48| ), and the v.e.v of the field A 
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plays the role of the physical mass m^ [0] . We can now write the equation for 
the spectrum. It can be obtained by writing small fluctuations around the 
master field a = ctq + -j^^<^y A = Aq + -^^A. The equations that one obtains 
are practically identical to ( p.24|) , apart a term Sa^x^X^ in coordinate space, 
that ensures 6axx = 0. It constrains the integral of the ip field to be zero. 
The eqs. of motion for the fluctuations are (see eq. (|3.24|) ) 



M2 



x{l — x) 



<^{x) 



— P I TTT T—rrz ^ ^{y) 



TV Jo {y — xY 
5X 



2x2(1 -X) 



2x(l-x) 



and 



dy^piv) = 



If we define 



F(m,r^) = P 



dx 



r"^ x{l — x) — w? 



(3.50) 



(3.51) 



(3.52) 



we get, for (7 = 0, the equation F{m,r^) = 0, that cannot be satisfied unless 
r^ = oo. If 5^ 7^ we can just use the constraint /q dxip{x) = to obtain 6X; 
using this equation back in the original equation we get 



M^ 



x(l — x) 



^{x) 



tp 



dy {x + y){2-x-y) 



+ 



with 



U^{M,r' 



-P 



-K Jo {y — xy 
x(l — x) 

dxdy 



2x2(1 -x) 



¥^{y) 

(3.53) 



{x + y){2-x-y) 



47r F(M,r2) 7o r2x(l-x)-M2 (y - xf 



viy) 

(3.54) 

This equation is very similar to the one in ( |3.24| ), where now the term 
with U^ takes the place of the last term with -^ in eq. ( |3.24| ). 
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4 QCD2 with quartic fermion interaction 

We can add to the Lagrangian of QCD2 a term quadratic in the fermion den- 
sity p without destroying the renormahzabihty of the theory. The request of 
Lorentz invariance leaves us only three terms: {ipipy , {ip'j^ipy and {ip'j f^ip)"^ . 
The resulting effective action will have a local term that is not linear in the 
p-densities. 

The model is described by the following Lagrangian: 



- /i(^^)'-/2(^75^)'-/3(^7m^)' 



(4.1) 



where the /« are three dimensionless coupling constants. In the light-cone 
gauge, we can eliminate the gauge field A°i to get a nonlocal term in the 
action 



S 



(fxd'^y \ ill) pip — rriQ ipip — ttiq ip'j^ip — fi{ipipy 



1 1 



- f2{^p75^Pr-f3ii^7,^r\s.y + ^9'r^r 



(4.2) 



We can write the effective action in the bilocal field formalism using as 
master field the (bilocal) matrix 



U 



PR P- 
P+ PL 



(4.3) 



We get (see [^ for the calculation of the Jacobian) 



2 
5e///iV = Tr{DU + i\ogU) + ^- j Sxd'^yG^yUllU, 



127-^12 

yx 



jd'xd^y[\[{ull)\[ull)' 



hUllUll + cUllUll 5.y 



(4.4) 



where 
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D{x,y) 



^xy 



(1) , (2) 



id- 



(1) (2) 



(4.5) 



V2 



and a = (/i + /2)A^ , h = 2f^N , c = (/i — /2)A^ are kept fixed as A^ — ^ cxd. 
The theory is chiral invariant if itlq = a = 0. 

We can move to Fourier space, and if we impose the saddle-point condition 

'' " (4.6) 



W\-q, -p) 







U=Uo 



we get 




D + jj-J {p,q)+g'6^'6'' 

(27r)2 (2^ 



d^k 1 



t/o^^(fc + p,-fc + g) 



(27r)2 F 
/ T^^-Sv^^^'^P'^ gi)(2vr)^5(^)(pi + gi - !> - g) 



with 



Mn 



af/o" + cf/|?2 



6^0^^ 



bU^' 



aU^'^ + cU^^ 



(4.7) 
(4.8) 



Notice that, iia = b + c = Mq is proportional to Uq^. 

Multiplying both sides of the saddle point equation with Uq {p, q) and 
choosing the usual translationally invariant ansatz Mo{p, q) = {2tt)'^6{p + q)Mo{p) 
we obtain 

{{D - A)Uof (p) + tS^' + g'6'' [ ^U','{p + k)U^\p) = (4.9) 



(2vr)^ 



where the matrix A is given by 

A'^ -- 



_fk_ 

(27r)2 



M^k) 



(4.10) 



If one compares eq. ( [4.9|) with eq.(13) of Ref. |2^ one can see that the only 
difference between the two is the shift: 



P-. ^ p^ — A^"^ 
p+ ^p+ - v42i 
moR ^ mon + A 



(4.11) 
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with moL, moR = (mg =F^o )/v^- The solution of the saddle point equation 
is given by 



-moL - A^'^ p 

1 

2p 



u,M = -2DAP') U , , . r(p '_) -,„„; 1 A- (4-12) 



with 



(4.13) 
r(p.)^^-^M (4.14) 
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We have changed the definition of r(p) with respect to Refs. [Q] and [^], 
so that we have replaced r(p_) <-^ — p_r(p_); the definition of the shifted 
momenta is p '_|_ = (p — a)± with a+ = A^^ , a„ = A^^ (see eq. ( [4.11| )). 



Now we have to enforce eq. ( |4.10| ) as consistency condition. A simple 
calculation shows that {A^'^ + A^^)/2 = 6m^^^/V2, {A'^^-A^'^)/2 = bmP-^ j ^[2, 
with 

5m(i) = {a + c) f -^ ^l"^' , (4.15) 

K m ( \ f ^^P 2i 1712 ,. ^„x 

dm} ' = [a — c) I - — — — - — (4.16) 

J (27r)2 2p^p^ — m^ + r(p_) + ie 

and we have defined the renormalized masses as rrii = rriQ + 5m^^\ and 
m^ = m^ — m2, so that we get the dependence of the coupling constants on 
the renormalization scale /Xq: 



27r 1^1 , f f4 



a + c 6m^^^ \m'^ — g'^ / Tc J 

'' - '"^ log M^ (4.17) 



a — c 6m^'^^ \im?—g'^/Ti^ 

A similar behaviour was already known in the context of the study of 
the CP'^^^ model with quarks, (see eq. (57) of [§], for instance), and as in 
eq. ( p.49| ) the only difference, due to the gauge kinetic term, is to shift the 
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squared mass of an amount 5'^/vt0. 

Furthermore, a careful analysis shows that after regularization A^"^ = A^^ = 0, 

due to the fact that U^^{p) is odd in p. 

In conclusion, the only effect on the master field of the interaction terms 
is to renormalize the mass, just as in the Bosonic theory and Uq is given by 



Uoip) 



-2i 



V2 



2p+p^ -m? + r(p_) +it\p+ + 2^r(p 



P- 

mi +9712 
V2 



(4.18) 



The equation for the fluctuation can be obtained as in |2^, where, using this 
formalism, we solved QCD2 with many flavors and with a chiral mass term 
m^^^ ipip + m^'^'' ijjj^ip . For the sake of simplicity we choose a = c and m^'^^ = 0, 
so that rriL = rriR = m/\^. Defining 



^/V2 



C^{p,-q) 






1 

m/V2 



(4.19) 

the form of the resulting equations suggests the following ansatz for the 
solutions 

¥P*^'(r,s_) 



"'t 



s\\C^ s_ 



Y,^--y)v^(^-) (4.20) 



so that finally we get after the usual rescaling of variables the correspondent 
of the 't Hooft equation for QCD2 with quartic fermion interaction: 



M2 



X 1 



X] 



9' 



vr Jo [y — xY 



dy 



+ 



a + c 



27r 



-771 — 



H m 

2-K 



vjy) 

X Jo I — y 
1 rl 



dy + 



1 



^{y) 



x(l — x) 



1 — X Jo y 

^ f{y) 



Jo y[l -y) 



v{y) 
dy] 

dy] 



(4.21) 



§ It is important to notice that in the case of a = c and m, 
scale of the spontaneous symmetry breaking 



(0 _ 




0, eq. (4.17) gives the 



6m — g /tt + /iQ exp 



27r 



We have a difference with respect to the analysis made in |l9] due to the gauge shift 
g'^/iT that appears in the correction Srri^ but that cancels out in the shifted mass 



M' 



dV^ 



g /tt giving substantially the same result. 



As in section 3, we can write the solution in the case b ^ 0, a + c ^ 0, g = 0. 
The previous equation becomes (we rescale r^ = rn?ii^) 



1 



/i 



27r 



x(l 
a + c /l 

\x Jo 1 
1 



if{x 



dy + 



fiy) 



b 

+ 2^ 



v{y)dy 



1 — X Jo y 

^ v{y) 



that gives 



if{x) 



x(l — x) Jo ^ ^"^ " ' Jo y{l — y) 
b A + pB + 2Bx{l 



dy\ 

dy 



X) 



27r 



li^ x{l 



X] 



(4.22) 



(4.23) 



where ^ = / 0(x), B = J (j){x)/x, p = {a + c)/b. This solution is remarkably 
similar to the one obtained by Fujita and Ogura [|T^] for the massive Thirring 
model |21| with Nf = 1 ; indeed, due to the equality 



{.l^.)ab{l'')cd = \d t 



ad ^bc 



(75) ad (75) fee 



(4.24) 



in the case of Nf = 1 our starting Lagrangian and the one of the Thirring 
model coincide for a = c = 6/2 = g*/4: where g* is the Thirring coupling 



constant (there is a different sign with respect to [|T8l due to the anticom- 
mutation of the fermions caused by T,aTl^jT^i and a factor of 2 that follows 
from our definition of the fermionic densities p ~ \/2^jj^jj). 
As in [jl8[ from eq. ( [4.23| ) we obtain the two equations 



A = — 



B = — 



b 

b 

2^ 



F{p'){A+[l 



/i 



4UU4b 



/i^ 



/i- 



^F{p')lA+\l + -\B) + {A + B)\oge 



/i- 



(4.25) 



where e is a cutoff needed to regularize the divergent equations, and 
F(/i2) = F(l,rVm2) (see JW^)- 

We must renormalize ( [4.25|) by adding a mass counterterm that verifies 



Sp^ (^(x) 



b_ 

2^ 



ei + 



£2 



x(l 



X] 



(4.26) 
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Adding (|4:22| ) to QOel ) we get 

b A + pB + e2 + {2B + ei)x{l 



(p{x) 



27r 



filx{l-x) -1 



(4.27) 



where /i|. = yU^ + (5/i^. In the following we will omit the suffix R for short. 

From ( |4.25| ) we observe that choosing ei = —b/27T ■ 2{A + B*) loge, 
£2 = -6/27r ■ p{A + B*) loge, with B* = B - b/27T ■ {A + 5*) loge we get 



^{x) 



b A + pB* + 2B* x(l - x) 



2tx p?' x{\ — x) — 1 
Now both A and B* are finite^ so we can solve the equations 



(4.28) 



A 
B* 



b 

2^ 
b p' 

2^^ 



f{p')Ia+ p + ^]b*] + ^b* 



p- 



p^ 



(4.29) 



■F{p') 



^+\P + 



B* 



getting 



b 

IT 



p^ , b 

P^ + ^ + ir 

4 Arc 



p^ 



F{p' 



1 



(4.30) 



Let us ffist examine the case p = 1. We can find bound states if ( [4.30| ) 
has a solution for < /i^ < 4. With the usual parametrization p = 2 sin6, 
< 9 < n/2, we have to solve the equation (A = —b/2n) 



A^ - A (3 - cos 26) + 
that, solved with respect to A, gives 



sin 29 
29 







(4.31) 



2A = 3 - cos 29 ± J (3 - cos 29^ - 4 



sin 29 
29 



(4.32) 



When 9 varies from to 7r/2, A„ goes from 1 to 0, whereas A+ goes from 4 
to 1; this gives the allowed range for the coupling constant in order to have 
bound states: < A < 4. In Fig. 3 we plot the two branches of 2A(^), and one 
can note the singularity at A = 1 of the function ^(A), that will be analyzed 
later. If A = 1 {b = — 27r) we have p^ = i.e. at the singularity we find 
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the well known massless pseudoscalar bound state of the massive Thirring 



model IjTj], that corresponds to the fundamental boson of the sine-Gordon 
theory |T^. 

Finally, we want to investigate the presence of tachyons in our model by 



searching a solution of ( [4.30|) in the range /i^ < 0. Putting /i^ = —Au"^ we get 



A2 _ 2A (1 - z/2) + /^^ I ""' - = (4.33) 



log(z/ + a/T 



u 



2} 



The equation has solutions iff A < (see Fig. 4); at A = 1 — z/q where 
uq ~ 1.59803 we have the 'lightest' tachyon, with mass /Xq = — 4z/q; for any 
value of A > eq. ( [4.33| ) doesn't have any solution, so we conclude that 
A > is the physical range of the coupling constant. In the tachyonic sector 
we have the following limits: 



1 



fi^{X) ^ -e"^ A ^ 0" 

/i2(A) ^ 2A A ^ -oo (4.34) 



For a generic p we have 



, 9 , r / ^NT sin 26* , , 

A^ - A [2 + p (1 - cos 29)] + — — = (4.35) 
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that, solved with respect to A, gives 



2A = 2 + p - p COS 20 ± i/(2 + p - p COS 2ef - 4 — -— (4.36) 

V 20 

Now the allowed range of A goes from 1+p— |l+p| to l + p+|l + p|: when 
< < 7r/2, A is always positive if p > —1, and always negative if p < —1. 
We can also see that for any p = implies A = 1 i.e. at the critical point 
A = 1 (6 = — 27r) we have a massless bound state. 

Furthermore, if p < pc = —1/3 the function 0{X) has a turning point at 
= ^0 > 0, as it can be seen from the expression of the argument of the 



square root in (|4.36| ) for small 



(2 + p - p cos 20f - 4 ^^^ ^ (2/3 + 2p) 40^ (4.37) 

20 

The behaviour of the argument of the square root in (|4.36|) is shown in Fig.5. 
We have forbidden values of if the function is negative. The value of the 
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turning point ^o increases for negative p and covers the whole range of 9 
at p = — 1 (Fig. 6); after that, the gap decreases in extension. The point 
p = — 1 is very important, because it is the supersymmetric point in the full 
theory with fermions coupled to bosons p[. At p = — 1 we can only have the 
massless bound state at the critical point A = 1, because ( [4.36|) has no real 
solutions except A = 1 (6' = 0). 

We can also see that near the value A = 1 the angle 9 has the behaviour 



3v^ 



;v/r^A A^r 



V3+yi 

:v/A^T A^r 



W^<! ^X f^-r . .. (4-38) 



if p = —1/3, and 



Vs-v^ 



"'^'^^^/a^^TsTsl^-^l ^^' <^-^^' 



if p > —1/3, that shows the transition at p = —1/3 and the singular be- 
haviour of ^(A) near A = 1. 

We can see more clearly what's happening if we investigate the tachyonic 
sector. The coupling constant is given by 



A = 1 - pz/2 ± 



z/Vl + z/2 



\ log(z/ + Vl + z^ ) 

In Fig. 7 we can see the behaviour of the argument of the square root in 
( [4.40|) : the function z/(A) has a turning point z/q > (for example vq ~ 1.59803 



if p = 1) that disappears at p < pc = —1/3, the critical point seen before. 

The special case p = contains no tachyons for any value of A, since it 
can be seen that (^4.40|) has no real solutions with nonzero v. 

Thus, it has turned out that in the tachyonic region A has always the op- 
posite sign of p: we have tachyons in our theory if pA < ; since A = — 6/27r, 
p = {a + c)/h we can conclude that the presence of tachyons depends only on 
a + c : we have tachions in the theory if a + c > 0. Our theory of self-coupled 
fermions is well defined provided that the coefficient of {ipipy in our starting 
lagrangian ( ^Tip is positive (i.e. /i < 0). 

In conclusion, in the free-tachyon sector we have a bound state for 
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2(1 + p) < A < if p < -1 and for < A < 2(1 + p) if p > 0. The mass of 
the bound state is given by M = 2mirsin^(A), with 9{\) given by (|4.36|) . If 
p = — 1 the only allowed value for the coupling constant is A = 1, for which 
the tachyon disappears and we only have a massless bound state; in the range 
—1 < p < we cannot have bound states without having also tachyons, so 
we won't examine this range of the parameter p (and the point pc = —1/3) 
far in detail. 

We finally observe that an interesting case is 6 = 0, a + c = 2/iA^ ^ 0, 
that corresponds to the Gross- Neveu model [|l9l if one identifies /i '^ —g^/2, 
g^ being the Gross-Neveu coupling constant (not to be confused with the 
gauge coupling), as used in ^^. One can see that the previous calculations 
lead to 

l = ^p^F(p) (4.41) 

that admits a bound state (note that such a bound state is never massless) 
if a + c < 0, given by the solution of (p = 2 sin 9) 

l = _^±f^tan0 (4.42) 
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If a + c > equation ( |4.41| ) has a solution in the tachyonic sector and the 



theory does not admit a ground state [|T9l ||20|. This confirms the results 
obtained before. 
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5 Complete gauge theory 

Finally, now that we've got experienced in using this formalism, we are ready 
to solve the full gauge theory (|2l] ): 



- 1^2(0^)^75^^ - Fi(0V)(V5^)' - F2(0V)(V575V')' - F3(0V)(^7mV')' 

(5.1) 

where we require the potentials to be polynomials in (j)^(f). Using the light 
cone gauge we are able to remove A'^ from the Lagrangian: 

c = ij{i^)ij + (f)\-d'^)(j)-v{(f)^(f))-Wi{(j)^(f))'ipij-W2{(j)^(j))'ip-f5'iP 

(5.2) 



- Fi(0V)(^^)2 _ F,{<p^){^p^,ijf - F3(0t0)(^^^^)2 + -g'r^^r 



This is our (nonlocal) starting Lagrangian, with J^ that is given by 

r = V'TVV^ + ^0^^" di (5.3) 

5.1 The 't Hooft equation for fermion-boson bound 
states 



We begin discussing the case 



W,{<P^<P)=moF (5.4) 

M^2(0V) = 

F,UU) = 



First, we eliminate the unnecessary degree of freedom tpi from the La- 
grangian (its kinetic term only contain (9_), and in the following we will 
define ■0 = 2iijjji (this is to remove boring \/2 terms out the Lagrangian); 



24 



then, to rewrite the whole interaction Lagrangian, we have to define four 
bilocal fields: 

p-y^Y.My)M^) (5.5) 

A 

a^y ^Y.^^A{y)M^) (5.6) 

A 

X^^ = E<^a(^)^^(^) (5-7) 

A 

^-y^Y.My)M^) (5.8) 

A 

In terms of the previous bilocal the interaction Lagrangian becomes: 

2 

Ant = ^^..(p^V"^ + (jy-i di I S a^y - x^"^ 5 r" + x'"^ S X"') (5.9) 

Now, it is a bit more difficult to get the Jacobian of the change of variables 

it is expressed by: 

J[p, cr, x,x] = J V2pV'ipV<P^Vcj) 6{p - ijij)S{a - 0V)'5(x - i'(p)S{x - 0^^) 

(5.10) 
By using auxiliary fields to exponentiate the (5's, we can calculate the Jaco- 
bian in the large-N limit: 

J = exp{-N STr log U) (5.11) 

with 

STr is the Supertrace, defined as 

STr 1^^ ^M = Tr(ml) - Tr(m2) (5.13) 

with rrii and /ij that are matrices with commuting and anticommuting ele- 
ments, respectively. This explains the choice of the — signs in ( [5.12| ). 
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By rescaling all fields by A^, keeping g'^N = g"^ constant we get the full action 
of the model: 



S/N = TrpgSTr(DU + ilogU) 

+ 2' J J2^^j2^^j2^^ ¥ U^{p„p,)W{k-p,,-k-p,) 

(5.14) 

with 
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D^Kp.q)=\^^ n^ 2 ° 2 (2^)''^^'^^ + ?) (5.15) 

\^ p'-m'oB ) 

"^^•(^' ^' ^) = ( (2p_ + A;.) (2p_ + M(2g^J- k_) ) ^^'^^^ 

The saddle point equation is obtained by imposing 

5S 



= (5.17) 

U=Uo 



5W'{-q, -p) 
that gives 

Ul\p,q) = U'M5^^^^\p + q) (5.1 



with 



rll/„N _ .„N -2^P- 



U^\P) = Poip) = , .2 .^ b-l ^ • (5-19) 

2p+p_ -M^-^^ + ze 



and (note the sign due to the Supertrace) 



- Ui\p) = aoip) = ^, |^_| . (5.20) 

The equation for the bound states can now be trivially obtained. We start 
from the equation for the fluctuations 

{l-26^yU^^{p,q) = -tg'U^^{p)Ul\~q) J ^^^^^ ^-^(-g^' -P' ^) 6U^^{k+p, -k+q) 

(5.21) 
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that is the generahzed equation for the 't Hooft "blob" , which has a leg made 
of a particle of type i and the other one made of a particle of type j. The 
6U^^ component leads to the 't Hooft equation 0, already obtained within 
this formalism in [^; the SU'^'^ component gives ( |3.19| ). Finally the other 
two components dU^"^ or SU"^^ give the equation for the scalar-fermion bound 
states: 



r'^ip{x) 



Ml Ml 



X 



(1-x) 



9' 



dy {x + y) 



^{x) - —P . 

7T Jo [y — xy Ix 



v{y) (5.22) 



recently written by K. Aoki 



5.2 Solving the most general Lagrangian 

Our method works also in the case of an arbitrary potential of the form 



v„. a' 









^1 "■' ^^ 

fii) (jn 
J n ^ 



ra=0 



n\ iV^+i 



(5.23) 
(5.24) 



Here the crucial difference is that we cannot remove anymore ipi as before 
and therefore we have to introduce a 3 x 3 bilocal field matrix: 



U 



( PR P- ~XR 
P+ PL -XL 

\ XL Xr -o- 



(5.25) 



As a consequence we get 
D{p,q) 



/ (1) , (2) 



V2 



\ 



P+ 





(1) (2) 





\ 



P^ - m^B ) 



{2nf6^^\p + q) (5.26) 



and the gauge part now couples the fermionic and the bosonic sectors. 
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We can also define 



V{Na) = N [mlB(T + V{a 

W,{Na) = m« + W,{a) 
A{a) = N[Fi{Na) + F2{Na)] 
B{a) = 2NF3{Na) 
C{a) = N[Fi{Na) - F2{Na)] 



(5.27) 
(5.28) 



and the L/R p-densities are defined (apart a y^ factor) as in [p3 |. 

If we impose tlie saddle point condition (|5.17|) to tlie resulting effective 
action with the usual translationally invariant ansatz we get (defining as in 



The self-interacting part can be written as 



(5.29) 



+ i^A{a) 



p'')" + (p''^" 



+ i3(a)pV'+C(a)p"p''(5.30) 



the gauge coupling part is given by the matrix W^ that has four nonzero 
elements: 



H^\U-p,q-k) 
H'\U;p,q;k) 
H'\U;p,q;k) 
H'%U-p,q-k) 

and the matrix Mq is given by 



{2p+k)U^^{k-q,-k-p) 

U^\k-q,-k-p) 

-{2q-k)U^^{k-q,-k-p) 

(2p + k) {2q - k)U^-\k -q,-k- p) 



Mn 






Cpf 



Bp]? 




\ 



Apf + Cpl 





-Ml' J 



(5.31) 



(5.32) 
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If Uq is translationally invariant it is easy to see that /Mq is a constant. 
It can be demonstrated that only the diagonal terms of the matrix A^^ are 
nonzero, so that we can absorb them in a mass renormalization. The master 
field Uq is block diagonal: 

P}^ Pi' 
Uo = I pf pf I (5.33) 



-ao 



For 2, j = 1, 2 we have 



(1) (2) 



™0 -^0 _ ^11 

p}^{p) = -2DM\ ^\ ^w^^r^) ,, I (5.34) 



with 



Df{p) 



2p+j9_ - (mf,^^ - mj)^^ + y2A")(m[,^^ + m[,^^ + ^2^22) + r(p_) + ie 

(5.35) 
and 



- f/f (p) ^ ao(p) = / 2 ^.33^^w ^^• (^-36) 

2p+P- - ("^05 + ^ ) + r(p-) + le 

The constants A*-' = /Mq"' are defined by using ( ^.32|) ; one can see that 
y4^2 = ^^21 _ g g^g before; the remaining constants must be regularized and 
then re-absorbed in a mass renormalization: 

V2A'' = 6m^'^ - 5m(2) 

V2A^^ = 6m^^^ + 5rrP'^ (5.37) 

2 
B 



^33 = 5m2 



so that, after the introduction of the renormalized masses m\ = rn^Q + 5m\ , 
nii = rriQ + ^m'-*-' , rnj^ = ml — m^ and of the shifted masses M^ = m^ — g'^/ir 
a = B,F we can write the propagators in a more familiar form: 
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Po{p) 



r," / mi -7112 
^^ ^/2 


P- 


2p+p_ - M2 - ^2|p_|/vrA + ze V p+ + 2^r(p_) 


mi+m,2 
V2 






'^^^^ 2p+p„-M|-^2|p_|/^^ + ,, 





(5.38) 
(5.39) 



It is now easy to get the equations for the fluctuations of the master field, 
by applying the same methods as in the previous section. Because of the 
presence of higher than quadratic terms in the densities appearing in the 
lagrangian, as in the 0{N) theory, we have to renormalize all the constants 
that couple two densities i.e. the scalar coupling V2, the Yukawa couplings 



w 



(1) 



,(2) 



and wi of (|5.24|) and the constants Uq , bo , Cq that did not need to 



be renormalized in the version of the model {QCD2 with quartic interaction 
term) studied in the previous section. This is our renormalization choice: 



1 ^-11^^ 



2" ^—^[V'FJ 
W1 + W2 _ W{{Ab) + W^{Ab) 

W1-W2 _ W{{AB)-W^iAB) 



V"{Ab) + ^A"{Ab) {a'pJ + (Af ) J + C"{AB)A'Mf 

+ A'iAB)A]} + C'iAB)Af (5.40) 
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+ A'{Ab)AJ,' + C'{Ab)A'p' 



V2 V2 

a = ^(Ab) b = B{Ab) c = C(Ab) 

Here the prime means derivative with respect to a, and we defined three 
(infinite) constants A^ = / o"o , A^ = J Pq. 

With the help of these definitions, we can write the form of the matrix 
6M, obtained by expanding the potential ( |5.3(]| ) up to second order in the 
fluctuations (compare with ( |5.32|) ): 



ri2 



6M 



\ 



y 5M33 j 

(5.41) 



30 



with 



5M' 



33_„,xrr33 , ^l+^2^^n^^l-^2^,,22 



v5U''' + 



-5U' 



(5.42) 



Our final result is the equations for the fluctuations in the most general 
case of a 2-dimensional gauge theory: 



l-28''Y8U'\p,q) = g'U',\p)Ul\-q) 



TJl 



dPk H'^{6U;p,q;k) 



k^. 



(fk 



U'o' (P) Ui' i-q) I T^6M^\k + p,-k + q) 



i2ny 



(5.43) 



These equations are the generalization of the ones already seen; the effect 
of the interacting potentials is to add inhomogeneous terms in the integral 
equations of motion. If one writes the entire set of equations, one can see 
that the independent degrees of freedom are the Fermion-Fermion and Boson- 
Boson fields p^^ = p_ = ipRipR: cr = 0V and the Fermion-Boson bound state 
Sxr = 'P^i'R- Since this last one is not subject to the interactions given by the 
additional potentials V , Wi and Fj, it does not develop an inhomogeneous 
term in his equations of motion. These are the full equations for the Fermion- 
Fermion and Boson-Boson bound states in the case of no chiral ipj^ip term 
rriR = rriL = itlf, Wi = w, W2 = (for the equation of the Fermion-Boson 
bound state, see ([5.22|)): 



Ml 



x(l — X 



ifpix) 



a + c 2/1 
-m ' 



2n 



+7^'mp 



27r 




n2 .1 



dy 



"^'^y^dy+ ' 



n Jo {y — xy 

^ VF{y) 



I — X Jo 



y 



^F{y)dy+ I -7^—jdy 



w niF 

271 x(l — x) Jo 



^Biy)dy 



y(l-y) 




(5.44) 
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and 



r^- 



M| 



x(l — x) 
1 fv 



9^ n /"^ '^?/ (a; + y)(2-2;-y) 



4x(l — x) \7T 



TT Jo {y — xy 4x(l — a;) 



VB{y) 



For f7 = we can give the solution of those equations: let us define 



a + c 

P = —r- ^ 



raph 



w wmp 

v = 



A= j\p{y)dy B=f'^^^dy C = j\B{y)dy 
Jo Jo y Jo 



so that 



ifpix) 



C + 7]B 



(5.45) 

(5.46) 
(5.47) 



Att r^ x{l — x) — m\ 

bml A + pB + aC + 2B x(l - x) 



(5.48) 



27r rpx{l — x)—m'p 

We then get a new set of equations that after the renormalization of B 



see (|4.25| )) are given by: 



A = — 

2tt 



B* 
C 



F{pb){C + ,,B-) 



A'Km'^Q 



(5.49) 



If 1 — v/4:T{m\F{pB) = then B* = 0; the only solution that doesn't imply 
also A = C = requires fi'^p = and 



A{1-X) = -^^^C 



27r 



(5.50) 



where we defined A = — 6/27r. 
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If 1 — v/4TTm'^^F{fiB) 7^ we can solve the matrix equations by obtaining 

A = 2fil{l- \)B* (5.51) 



sin 26b 
29r 



+ 



Airm 



B . 



\'^-\[2 + p{1-cos29f)] 



sin 26f 



29, 



C 



w 



wrriF 
47rm| 



B* 



(5.52) 



2™,,:^-"^^''-HiFl = °(^^^^' 



C 



where we made the usual parametrization /x^ = 2 sin 9a, < 9a < 7r/2, 
a =F,B. 

The equations depend on two constants, that we can fix by normalizing 
ipp and ips] requiring A = C = 1 we get 



1 



+ 



wmp Up 



4:TTm\ Sirm'^^ 1 — A 



F(^ 



B 



1 = -2\F{fXF) 



Hf w 

4 Attitlf 



1 



1-A 

2 



(5.54) 
(5.55) 



The second equation can be solved as a function of A 



2A 



2T[rai 



+ p(l-cos 29f) 



± 



\ 



w 



2nmi 



p (1 - cos 29f) - 4 



sin 29f 
29f 



w 



2Trmi 



(5.56) 



We will not proceed further in solving those equations, but we can give a 
perturbative estimate of the shift in 9a {w small): 
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wniF 1 1 — cos 29f 

2S9p = ^ _ 1-^ 

2iTmF h{9F) — Xp sin 29f 

where h{x) = d/dx{sm 2x/2a;), 9a = 9a + 59a. 



(5.57) 
(5.58) 
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6 Conclusions 

We have shown the vahdity and the flexibihty of the bilocal field method re- 
cently proposed ||23|; we have reproduced via the path-integral approach all 
the known results about models as 0{N) and massive Gross- Neveu model in 
the large- A^ limit, and we have generalized the 't Hooft equations for QCD2 
and scalar QCD2 to the case of self-interacting fields. 

Our formalism is a kind of generalization of the one used in the context 
of the 0{N) models [§], but also applies to all cases in which the field is a 
'vector' in the colour index, independently if the vector index belongs to an 
internal symmetry group (such as 0{N)) or to a gauge group. 

We could also conclude that the two-dimensional gauge models with mat- 
ter in the fundamental representation are as trivial, from the point of view 
of the large-N expansion, as the 0{N) and CP^~^ models, as we just have 
to pay the price of introducing a bilocal colour-singlet field instead of a local 
one. This is not the case of matter in the adjoint representation, in which 
one can construct an infinite number of multilocal fields that all contribute 
for large A^. Such a structure enlarges the difficulties but gives a richer spec- 



trum of states, both in the case of QCD2 with fermionic [^ or bosonic p5 
matter in adjoint representation. 
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Fig.l Graphic resolution of the equation x=a*Sin(x) for a>l, a<l 




Fig. 2 Graphic resolution of the equation x = a*f (x) for a>l, a<l 
f(x) = (l+x''2) '^-1/2 * Log(x+sqrt (l+x''2) ) 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9401086vl 




Fig. 3 Coupling constant as a function of the mass angle x (eq. (4.32) 
In the bound state region M=2m*sin (x/2) 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9401086vl 




Fig. 4 Coupling constant as a function of the z variable (see eq. (4.33) 

In the tachyon region M'^2=-4m'^2*z'^2 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9401086vl 




Fig. 5 Plot of the argument of the square root in eq. (4.36) 

for rho = 1,-1/3,-3/5,-1,-5 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9401086vl 




Fig. 6 Plot of the argument of the square root in eq. (4.36) for rho = -1 



This figure "figl-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9401086vl 




Fig. 7 Plot of the argument of the square root in eq. (4.40) 

for rho = -1,-1/3, 0, 1 



This figure "figl-7.png" is available in "png" format from: 
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